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, Abstract 



We show a surprising connection between known integrable Hamil- 
tonian systems with quartic potential and the stationary flows of some 
I coupled KdV systems related to fourth order Lax operators. In partic- 

ular, we present a connection between the Hirota-Satsuma coupled KdV 
system and (a generalisation of) the 1:6:1 integrable case quartic poten- 
tial. A generalisation of the 1:6:8 case is similarly related to a different 
Oh! (but gauge related) fourth order Lax operator. 

We exploit this connection to derive a Lax representation for each of 
these integrable systems. In this context a canonical transformation is 
derived through a gauge transformation. 



X. 

■ 1 Introduction 



In a surprising connection was shown between integrable cases of the Henon- 
Heiles system and the stationary flows of some known integrable PDEs : the 
5th order KdV, Sawada-Kotera and Kaup-Kupershmidt equations, which are 
integrable through second and third order differential Lax operators. This gave 
rise to a matrix spectral problem for each of the integrable Henon-Heiles systems. 
Thus, integrable cubic potentials are associated with second and third order 
matrix Lax operators. 

In Q the starting point was the Henon-Heiles system, using some elementary, 
ad hoc calculations. An alternative approach is to start with the stationary flow 
of a known integrable PDE and to find the appropriate Henon-Heiles coordi- 
nates. The most natural way of doing this is through the Hamiltonian structure 
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of the PDE. As motivation, this is presented for the Henon-Heiles system in 
section 2. 

In this note we exhibit a similar connection between some stationary flows 
associated with fourth order Lax operators and generalisations of some inte- 
grable Hamiltonian systems with quartic potentials: 

H=^{pi+ pI) + aqf + hqlql + cqi (1) 

There are 4 nontrivial cases which are integrable: 

1. a:&:c=l:2:l, 

2. a : & : c = 1 : 12 : 16 , 

3. a:fe:c=l:6:l, 

4. a:&:c=l:6:8. 

Various inverse square terms can be added without destroying complete integra- 
bility. Cases (|l|), (^) and (||) are separable in respectively polar, parabolic and 
Cartesian coordinates. Case is one of the higher degree polynomial exam- 
ples given in |^ . Case (Q) is much more complicated and has only recently been 
shown to be separable (in a generalised sense) [3|, using a Painleve expansion. 

Our main result is to relate cases (||) and (4) to some coupled KdV equa- 
tions associated with fourth order Lax operators, which are themselves related 
through a Miura map. The connection with integrable PDEs gives a Lax rep- 
resentation for these finite dimensional systems, together with similarity and 
canonical transformation, thus giving an explanation of the separation coordi- 
nates found in ^ . 



2 Stationary Flows and Integrable Henon-Heiles 
Systems 

It was shown in [Q that the general (not necessarily integrable) Henon-Heiles 
system can be related to the stationary flow of: 

ut = {d^ + 8aud + Aaua,) Su J^^^ " ^"^) ■ (2) 

For this stationary flow, the gradient of the above Hamiltonian is in the kernel 
of the third order Hamiltonian structure. Thus we may write: 




= A, where [d^ + 8aud + iau^) A^O. 
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We now set A — ay^ to find: 

yiVxx + '2auy)x + iyxiyxx + 2auy) = 0. 

Setting yxx + 2auy = 7, we solve the first order equation for 7 to get 7 — 2ky~^, 
where fc is a constant. We now have: 

yxx + 2auy = 2ky~^ , 
which, for a = —a, are Lagrangian, with: 

^^=\{ul + yl) + \bu^ - auy"^ - ky'"^ . 

The standard Legendre transformation now renders a natural Hamiltonian sys- 
tem, which is just the usual generalisation of the Henon-Heiles system: 

qi=u, q2^y, P2 = yx, 

h^^{pl+ pI) - ^bqf + aqiql + kq-^ . 

Thus the Hamiltonian structure of (^) gave us a natural way of defining some 
interesting coordinates, giving us the Henon-Heiles representation of the station- 
ary flow. For the integrable cases it is possible to use the Lax representation of 

to derive a matrix Lax representation for the corresponding Henon-Heiles 
system Q, thus proving the complete integrability of the latter. The gauge 
equivalence of the Sawada-Kotera and Kaup-Kupershmidt equations leads to 
a canonical transformation between the corresponding Henon-Heiles systems 
(a/6 =-1/6 and a/6 =-1/16). 

In this paper, we obtain similar results for some quartic potentials. 

3 Fourth Order Operators 

We start with the self adjoint fourth order operator, which we write in factorised 
form: 

L^^{d + V,){d + V2){d~V2){d-Vi). (3) 

This can be written as the product of two second order operators: 

L,^{d^ + fd + fx+g){d^-fd + g), 
to give the Miura map: 

/ = t;i+W2, 9^viV2~vix. (4) 
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The Lax equation: 

Lit = [Afi,Li] 

where: 

Ah = 2d' + ^i2g - - f)d + ^{2g, ~ - 2//,), 
gives the following coupled KdV system: 

ft = -i(2/,,, + 3//,, + 3/J-3/2/. + 6/<?, + 65/,), 
gt = j{2gxxx + l2ggx + 6fgx^ + 12gfxx + 18fxgx-6pgx 

~^^fxxxxx ~t" ^ffxxx ~t~ ^8 fx f XX 6/ /xx ^ffx^' 



(5) 



Remark 1 T/ie resulting coupled KdV system is simpler in co-ordinates r and 
s, corresponding to the operator: 

Lsym = 9"* + drd + s, 

but the Miura map and Hamiltonian operator are more complicated, so not well 
suited to our purposes. 

A rotation of factors in leads to another operator: 

L2^{d-Vi){d + Vi){d + V2){d-V2), 

which can be written as: 

L2^{d^ + u + ip){d^ +u-if>), (6) 

with the Miura map: 

U =^{vix - V2x - vl - vl), (p^]^{vix+V2x-vl+vl). (7) 

The operator (^) was foimd in Q to be the Lax operator for the Hirota-Satsuma 
system: 

Ut = -Uxxx + iuUx-Q^P^Px, (fit = -Vxxx -^U(px, (8) 

which corresponds to the time evolution operator: 

M2 = 2d' + 3ud+^iux-2^x)- 

The isospectral flows of the general fourth order operator are bi-hamiltonian, 
but the first Hamiltonian structure does not reduce to the 2— component sub- 
manifolds corresponding to our two operators. As in the case of our third order 
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operators, both the /, g and u, if hierarchies are related to the same modified 
hierarchy for the functions vi and W2, which has Hamiltonian form: 



The Miura maps (^,0) transport this structure respectively onto: 

-2d -d^ -fd-u 



B/ 
and: 



52 -fd + (2/, + 2g- p)d + + 5x - fU) 



1 a3 



d'^ + 2ud + Ux 2tpd + tpx 



2ipd + ipx ^d^ + 2ud + Ux 



These are reductions of the general second Hamiltonian structure associated 
with the general fourth order Lax operator. With these, equations (||) and 
(||) are respectively generated hy Q = —\{fx + + '^fdx — — 4(^2) and 
7i = ^u^ — (f"^, both of which are pulled back to the modified Hamiltonian 

hrnod = -| {vfx + ^2x + '^'l + "^1) " | {'V\xV2x + v\vix - v\v2x " W^V^) , which 

evidently generates a coupled MKdV system. 

3.1 Stationary Flows and Related Quartic Potentials 

We now consider the stationary flows: 

B/JfS = and B„(5uH = 0, 

which we show to be respectively related to generalisations of cases (|) and (|) 
of the quartic potentials (|l|). 

3.1.1 The 1:6:8 Potential 

The stationary flow in the (/, g) space is given by: 

{2A + + !B)x = (9) 
Axx - I Ax + Bxxx + {2fx + 2g - f)Bx + {fxx + gx - ffx)B = 0, (10) 

where, for the Hamiltonian Q, 

A = 5fg=^{fxx-2gx-2f + Afg), (11) 

B = 5gg = ^ifx + f + 2g). (12) 

Equation (||) gives A — ^{K — fB — Bx), where JsT is a constant of integration, 
and: 



5 



B^^^ + 2(/, + 2g~ + + 2g^ - ff,)B = 0, 



which has solution B = ay^, where 



y.. + ilfx - \f + g)y ^ Ly-\ (13) 



Using (^2|) we obtain a formula for 



which, from (11), ( p^ ) and the definition of K, gives: 

- 2f + iay^f = K, 
3 

ya:x + (ay^ - -^f)y - Ly-"^- 
When a = — i the above equations are Lagrangian with: 

Cu,a) = \ul + vl) + Yy + e/'y' + 8/') + Kf ^y- 

With canonical co-ordinates: 

Qi = y, Q2 = /, Pi = yx, P2 = fx, 



this gives the Hamiltonian: 

= ^(A' + Pi) -Y^iQt + ^QlQl + - KQ2 + ^Q^\ (14) 

which is one of the integrable generalisations of case of (|l|) |^ . 



3.1.2 The 1:6:1 Potential 

The stationary flow in the {u, ip) space is given by: 

1 r, 

{-d^ + 2ud + u^)A+{2(pd + (p^)B = 0, 
1 Q 

{2(pd + (p^)A+{-d^ + 2ud + u^)B = 0, 

where, for the simple Hamiltonian Ti, we have A — Su'H = u and B = S^Ti 
■—2ip. It is straightforward to find quadratic forms for A and B: 

A^ai^^+X^), i? = a(x'-^'), 

where: 



Xxx + iu + ip)x = kx 



-3 
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Remark 2 The motivation for this representation of A and B is found in the 



squared eigenfunction representation of L2O 
second order differential system. 

Upon substituting u and (p we get: 

Xxx + ^0(31/"^ +x^)x 



— zO, when written as a 2 x 2, 



which is Lagrangian with: 



= \{€ + Xl) " + + X') + kx~^). 



With the canonical co-ordinates: 

qi pi = Vx, 92 = X, P2 

this gives (with a — I) the Hamiltonian: 

which is one of the integrable generalisations of case (||) of (|l|) |^ 



lipl +pI) + liqt + ^qlql + qt) + \{iqi^ + kq^^), 



(15) 



3.2 The Lax Representations 

The Lax representations Lu = [Mi, Li] for the PDEs can be re- written in zero 
curvature form: 



where: 



( 1 

-.9 / 



V z 





1 
1 

-.9 - !x -f ) 



U2 



( 1 

ip — u 



\ z Q 



\ 

1 
1 

—If — u J 



and V^i given by more complicated formulae. The stationary flows arc given by 
the Lax representations Vix = [Ui, Vi]. which can be written in terms of qi,Pi 
and Qi, Pi. Whilst Ui are very simple, we need to use the equations of motion 
generated by /i(/.g) and h(^u,ip) to eliminate second and higher derivatives in Vi. 
When written in this way the Vi are the Lax matrices for the equations of motion 
and can be used to generate the constants of motion. They are given by: 



Vi 



4 
021 

2Q2Z 



Q2QI 



- \QiPi 

2z 




2O2 

033 
043 



Ql 



2 


-\Ql 

0-44 



7 



v. 



where: 



2z - 2p{ - 2lq^^ 2qipi 




-2qipi 2ql 

2z ~2q2P2 ^H2 

2z- 2pl - 2kq^^ 2q2P2 ) 



2gi 



021 = 

033 = 

043 = 

044 = 



1 

16 



Q1Q2 



L 

4q! 



J^i + \Q1Q2P1 + 2z, 



1 

16 



Q1Q2 



L 
4Q! 



1 



1 



:QlQ2Pl + '2z, 



lQ2Qi - \qiPi + \k. 



The first integrals are given by the characteristic equations: 

det(yi - ^il) = /i^ + i (L - AK^)n^ + det Vi , 
8 

det(V2-Ai^) = // + 4(A: + 0M^ + detT^2, 



(16) 
(17) 



where 



-16z^ - 8/1 



1 



1 



det(Vi) 

Aei{V2) = -16z^ + 32/i(„_^)Z^ - 16fc(„,<^)Z + 16Zfc. 
Here and are given by ( p^p^ and: 



2\2 



Hf,a) 



Pi + ^{QlQl + QiQt) - -AQ\Pi + QtPi) + ^Ql + P1P2Q2QI - ^QlQlPf 

16 4 64 2 



-\{Q\Q2K + LQl) - 2{QlK^ - ^) 
-miPiP2K - Q2P?K - ^^). 



1 



L2 



QiQiK) 



P1P2 + 2(9192 + 91(72) 



Qi 

kql 



Ik 

mi 



-^+lql + 



The Lax equations for Vi are known to Unearise on the Jacobi variety of the 
algebraic curves defined by (|l6|) and ( |l7|) , which gives one of the standard inte- 
gration procedures (see, for instance, Theorem 1, page 67 of [^) 

3.3 Gauge and Canonical Transformations 

The matrices Ui, Vi (in the PDE variables) are related by gauge transformation: 



U2 = AUiA-^ + A.^A- 



V2 = AViA-^+AtA- 
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with A given by: 



-Vi 


1 







-VIV2 




1 













1 


z 





-V1V2 





v 

whenever the Miura maps and (j^) are satisfied. In terms of and their 
derivatives, the Miura maps can be rearranged to give: 



2,2 
ll + 

2 I 1 2 

92 + g^l 



2g2za; 
3(72 

2gl3;z 

3(7i 



^^92 "^^ 
2, -4 



We can then write vi and V2 in terms of qi and (72- 



q2x a qix 

vi — 1 — 5- and V2 — 

92 92 91 



91 



where = —k and 6^ = — L We can similarly write Vi in terms of Qi and Q2'- 



1^ I 'Sia: 
= 77(^2 + 



1 c 



"q7 



2Qf' 



where c = 8(6 — a). 

We use these formulae, together with the Miura maps and the defini- 
tions of the canonical coordinates, to explicitly write down the canonical trans- 
formation: 



Qi = 2T5, 

Q2 = ^-^ 
91 92 

and its inverse: 

1 



9i 



91 

92 



-r 



1/2 



1, 



,1/2 



9? 



Pi 



Pl _ P2 

91 



a 

92 ' ql 
2 



ql 



-^1/2 _ £^-1/2 

4 



D 2 2 , pI Pi 2ap2 26pi 
-f2 — 9i - 92 + — 

92 



9? 



9| 



9? 



Pi c 
^^1 = 2 1 ^"^2 - — 



1 fl^ 
1 /I 



J52 



2 V2 



:02 



Oi 2Q2 
^ Oi ^ 2Q2 



ri^^ - 26r 



-1/2 



U/2 



2aT 



-1/2 



9| 



4 ' 
91 



where 



T 

r± 



2 I 2 
9i +92 



9192 



P1P2 



±2P2 + Q^+2«f-^T 



&P2 api o-b 
91 92 9192 
4F2 ^Q^p^ 2cQ2 4cPi 



c 
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The original parameters k, I, K, L are related to a and b by: 

k = -a^, I = -52, L = -16(6 - a)2, K = -2(a + 6). 



Remark 3 When a ^ b = 0. the case 1:6:1 separates in Cartesian co- 
ordinates and the canonical transformation between cases 1:6:1 and 1:6:8, 
constructed in this paper, gives the separation co-ordinates for 1:6:8 obtained 
in by a Painleve approach. Some sort of Lax representation is given in 
, but this involves the square root of k(^f gj and its derivatives, so does not 
give a clear derivation of the equations of motion as a consequence of the Lax 
equations. 

4 Conclusions 

The whole essence of our approach is to give a systematic connection between 
completely integrable, finite dimensional Hamiltonian systems and integrable 
(in the soliton theory sense) PDEs. This result is, itself, interesting but, more 
importantly, gives a systematic construction of a matrix Lax pair for the finite 
dimensional system and (as a bonus) gives a straightforward construction of 
canonical transformations (via gauge transformations in the PDE framework). 

We just presented the most interesting cases ((||) and (||)) of (generalisations 
of) Hamiltonian (|l]) in this letter, but the method is quite general. Case (^ 
can be obtained from the KdV hierarchy, since it is just a special case of the 
Garnier system. Case (|]) is a special case of a class studied in associated 
with energy dependent Schrodinger operators. Case (|^) can also be obtained by 
considering a fifth order Lax operator 

A more complete treatment will be presented elsewhere. 
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